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Atomic supersymmetry is a quantum-mcchanical supersymmetry connect- 
ing the properties of different atoms and ions. In this talk, I provide a short 
description of some established results in the subject and discuss a few re- 
cent developments, including the extension to parabolic coordinates and the 
calculation of Stark maps using supersymmetry-based models. 


1. Introduction 

It can often be impractical to find exact solutions to the equations determining 
the properties of a physical system. The identification and use of symmetries is 
one method that can be useful in the search for a mathematically simpler but 
physically sufficient description. 

An example is provided by the behavior of the valence electron in alkali-metal 
atoms, which is governed in detail by the solution to an involved many-body prob- 
lem. The essential physics of this situation can largely be contained in a single- 
electron model with an effective central potential [1]. This talk concerns symmetiy 
issues involved in this approach, in particular, the role played by supersymmetry 
in a realistic central-potential approximation. 

The first part of this talk (sections 2-5) provides some background and an 
overview of established results in atomic supersymmetry using spherical coordi- 
nates. The second part (sections 6-8) discusses recently developed extensions of 
these ideas to other coordinate systems and to the supersymmetry-based descrip- 
tion of the Stark effect. More information on the approach used, its relation to 
other methods, and recent developments can be found in the references cited. 


2. Supersymmetric Quantum Mechanics 

A quantum-mechanical system is called supersymmetric [2] if its hamiltonian 
H commutes with N supersymmetry generators Q } satisfying {Qj,Qk} = &jkH. 
The superalgebra generated by H and Q j is denoted sqm(N). Of interest here is the 
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special case sqm(2), for which it is convenient to consider the linear combination 
Q = ( Qi +iQ2)/V2 satisfying H = {Q,Q^}. 

A simple representation of sqm(2) sets 

«“(a o) - *=(t L) • < 2 -» 

where A = —id x — iU ' / 2 with U' = dU/dx for some function U = U(x). The 
components satisfy 


J2 

■^±^±n = [ — = f n^in ? (2-2) 

where V±{x) = U' 2 / 4 + U" /2. 

If the supersymmetry is unbroken, the energy of the ground state is zero. The 
eigenspectra of H± are degenerate except for this state, which is associated with 
H + . The spectra of H+ and H— are called, respectively, the bosonic and fermionic 
stacks. Degenerate states in the bosonic and fermionic stacks are mapped into one 
another by the supersymmetry generators Q , QK 


3. Atomic Supersymmetry 


The Schrodinger equation for the hydrogen atom separates in spherical polar 
coordinates. The radial piece can be written in atomic units as 


r d 2 1 1(1+1) 1 

[ dy 2 - 2 0 En]Xnl(y) — 0 , 


r 


(3.1) 


where y = 2 r, E n = —1/2 n 2 , and x n j(2r) = rR n i(r). The radial wave functions 
are 

2 


Rnl(r) = — 


rr 


T(r 




where the Sonine-Laguerre polynomials are defined by 

T(n + o + l) 


lW(x) = Vf -iV 

" u p! ( n — p)- T(p + a + 1) 


(3.3) 


The idea is to identify the radial equation (3.1) for fixed l with the hamiltonian 
H+ of a supersymmetric quantum mechanics, as in Eq. (2.2). One can then 
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construct the supersymmetry generator Q, determine H-, and seek a physical 
interpretation [3]. 

The supersymmetry is determined by specifying the function U. For the 
Coulomb problem, it is U(y) = y/(l + 1) — 2(1 + l)lny. The partner hamiltonian 
H- can be shown to have the same form as H+ except that l is replaced with 
(/ + 1) in the 1 jy 2 term. This means that the eigenfunctions of H _ at fixed l are 
R n ,l+i w hh n >2. These form a complete and orthonormal set. 

To understand the physical interpretation, recall that H _ has the same spec- 
trum as H + , excluding the lowest state. This suggests H_ describes a physical 
system that looks hydrogenic except that the ground state is inaccessible. The 
idea is to attribute this to the Pauli principle. The partner system is then an 
atom with the lowest state filled with electrons, and H _ is an effective hamilto- 
nian for the valence electron. For example, if 1 = 0 the spectrum of H + spans the 
s orbitals of hydrogen and the spectrum of H— represents the s orbitals accessible 
to the valence electron in lithium. Evidently, the supersymmetry holds only if the 
non-hydrogenic electron interactions are disregarded. The incorporation of these 
symmetry-breaking effects is described in section 5. 

In the exact-symmetry limit, the supersymmetry can be extended. For ex- 
ample, the s orbitals of lithium and sodium can also be viewed as supersymmetric 
partners. There are supersymmetric connections of this sort among atoms and 
ions across the periodic table, all physically incorporating the effects of the Pauli 
principle and mathematically implemented by integer shifts in the angular quan- 
tum number l. 

Ref. [3] presents more details of the construction and discusses the experi- 
mental support for this atomic supersymmetry. 


4. Oscillator Formulation 

The oscillator formulation of atomic supersymmetry outlined in section 3 
is closely related to harmonic oscillators. Indeed, connections exist in arbitrary 
dimensions among the radial equations for the Coulomb problem, the harmonic 
oscillator, and the two sqm(2) supersymmetric partners of these systems [4]. 

The radial equation for the d-dimensional Coulomb problem can be written 
in atomic units as 

[~ - t + (f + 7>( | )2 + 7 + 1} - lg,„] v M (y) = 0 , (4.1) 

where y = 2 r, E dn = -l/ 2 (n -f q) 2 , 7 = (d - 3)/2, and the radial wave functions 
are 

Vd nl(y) = c dn i y f+7+1 exp(— y/2(n + 'y))L i ^^ 1 1 + 1) (y/(n + 7 )) , (4.2) 


445 


with Cdnl a normalization constant. If the operator in Eq. (4.1) is interpreted as 
H+ in an sqm(2) supersymmetry, the partner H- can be shown to have the same 
form but with / replaced by 1 + 1. 

The radial equation for the D-dimensional harmonic oscillator with unit an- 
gular frequency in atomic units is 

- y + + r)( y 2 + r + 1} - E dw ] V dnl (Y) = 0 , (4.3) 

where Y is the radial variable, Edn — 2iV + 2T + 3, T = (D — 3)/2, and the radial 
wave functions are 

v DNL (Y ) = c DNL F L+r+1 exp(-r 2 /2)i ( N L /t-t}2 2) (^ 2 ) » ( 4 - 4 ) 

where Cdnl is a normalization constant. The sqm(2) partner fZ_ to the operator 
in Eq. (4.3) again has the same form but with L replaced by L + 1. 

A state in the d- dimensional Coulomb problem can be related by a one- 
parameter mapping to a state in a D- dimensional oscillator [4]. The eigenfunctions 
v dnl and Vdnl are related by Urf„/((n+ 7 )F 2 ) oc Y 1/2 Vdnl(Y ) with N = 2n-2+A, 
L = 2/ + A, and even D = 2d — 2 — 2A. The parameter A is an integer. Restrictions 
on the possible choices of the parameter, the dimensionality, and the quantum 
numbers can arise if several Coulomb-problem states from one system are required 
to be mapped to a specified harmonic oscillator. The map can be combined with 
the supersymmetric maps shifting the angular quantum numbers, resulting in a 
commutative diagram between states of the four systems. The reader is referred 
to [4] for further details. 


5. Supersymmetry-Based Quantum-Defect Theory 

The energy eigenvalues of the valence electron in an alkali-metal atom are 
shifted relative to the hydrogenic values by interactions with the core. The 
eigenvalues can be expressed as the Rydberg series [5] E n > = —1/2 n* 2 with 
n* = n - 6(n,l ), where for given l the quantum defect 8(n, 1) rapidly approaches 
an asymptotic value 6(1) as n increases. These shifts break the supersymmetries 
described in section 3. 

Despite this, it is possible to develop a central-potential model incorporat- 
ing supersymmetric features that has analytical wavefunctions as solutions and 
eigenenergies given by the Rydberg series [6]. The theory is defined by a radial 
equation obtained from Eq. (3.1) by replacing n,l,E n with n* ,1* , E n * . Here, l* 
is a new angular quantum number given by /* = / — 8(1) + 1(1), where 1(1) is 
an integer shift characteristic of supersymmetry. The eigenfunctions are given by 
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the replacement of n and l in Eq. (3.2) with n* and /*. The Sonine-Laguerre 
polynomials still appear because n* — l* — 1 = n — l — 1(1) — l remains an integer. 
For asymptotic quantum defects <$(/), these eigenfunctions form an orthonormal 
and complete set. 

Ref. [6] contains more details about this model and demonstrates that the 
analytical eigenfunctions yield transition probabilities that are in agreement with 
experiment and with accepted values [7], Related developments, other than those 
to be discussed in the remainder of this paper, include the following. 

• The eigenfunctions of the model have been used as trial wavefunctions for 
detailed atomic calculations [8], 

• The validity of the theory at short distances has been examined by investi- 
gating predictions for the fine structure of alkali-metal atoms [9], Certain key 
features observed experimentally are reproduced and the Lande semiempirical 
formula naturally appears in the model. 

• The model has been used to examine transition probabilities for other atoms, 
notably alkaline-earth ions [10]. 

• The mathematical structure of the model has been shown to be connected 
with parastatistics [11]. 


6. Separable-Coordinate Supersymmetries 

The Schrodinger equation with a Coulomb potential V(r ) = — 1/r separates 
directly (no modulation factor) in four different coordinate systems: spherical, 
conical, prolate-spheroidal, and parabolic. Each separation results in three ordi- 
nary differential equations, so a priori there are twelve candidate equations for the 
role of a supersymmetric partner hamiltonian in an sqm(2) realization. However, 
additional constraints arise from the structure, of Eq. (2.2): to play the role of 

an equation must be expressed in the appropriate form and the eigenspec- 
trum must involve a tower of states labeled by the appropriate eigenvalue. Some 
results of an examination of possible sqm(2) supersymmetries involving the twelve 
candidate equations are described in this section. More details are provided in 
ref. [12]. 

First, consider spherical coordinates. The eigenfunctions for the Coulomb 
problem separate into a product of radial wave functions given in Eq. (3.2) and 
the spherical harmonics, which in turn decompose into products of associated Leg- 
endre functions and exponentials. As discussed in section 3, the radial equation 
does admit a supersymmetric partner. The equation for the associated Legendre 
functions can be rewritten in the form of Eq. (2.2) by admitting a modulation 
factor, but the ensuing equation lacks an eigenvalue representing an infinite tower 
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of states. The remaining azimuthal equation is uninteresting because it has zero 
potential. Therefore, the only interesting sqm(2) supersymmetry in spherical co- 
ordinates appears in the radial equation. 

In conical coordinates, the coordinate surfaces are spheres and cones of ellip- 
tic cross section along the 2 and x axes. One coordinate is just the radial variable, 
and the eigenfunctions for hydrogen are products of Lame functions with the same 
radial wavefunctions as in spherical coordinates. The radial equation has the su- 
persymmetry described in section 3 and the same physical interpretation applies, 
although the complete eigenfunctions for the supersymmetric partner states in- 
volve Lame functions rather than the spherical harmonics. Each of the other two 
separated equations can be rewritten without a linear derivative term as in Eq. 
(2.2), but no interesting supersymmetry exists. 

The Schrodinger equation for the Coulomb problem also separates in prolate- 
spheroidal coordinates. It can be shown that none of the three resulting equations 
admits an interesting sqm(2) structure. 

Finally, consider parabolic coordinates. These are defined by p — r + z, 
<7 = r — z, <j> = tan _1 (y/x), where r is the radial coordinate. The hydrogen eigen- 
functions in parabolic coordinates can be written as the product u(p)v(cr)<f>(<fi). 
The equation determining <&(</>) is uninteresting in the present context. The equa- 
tion for u(p) can be expressed in the form of an H+ by extracting a modulation 
factor u(p) = p -1 / 2 Xi(p). The ensuing equation for xi is 


r d 2 1 — m 2 
dp 2 4 p 2 


+ f + f]* ( ') = °’ 

p & 


( 6 . 1 ) 


where E is the energy and m is the magnetic quantum number. A similar analysis 
can be carried out for u(cr), resulting in an identical equation but with Xi 5 P , and 
Z\ replaced by X 2 , cr, and Z 2 . The separation constants Z\ and Z 2 are required 
to satisfy Z\ + Z 2 = 1. The eigensolutions have the form 


u 


n\n?m 


(p) 




f+ i \- 


T(n 1 + 1) 

T(ni + m + 1) 


— L 

e ^ 


( 6 . 2 ) 


with a similar expression for v„ in2m (a). Here, ni = Z\ je — (m + l)/2, n 2 = 
Z2/e — (m + l)/2 are nonnegative integers, and e = \J—2E for E n — —1/2 n 2 with 
n = rc 1 + ri2 + m + 1 . 

The equations for xi, X 2 admit supersymmetric partners. The functions U(p), 
U(<t) determining the supersymmetry are U(p) = 2Zip/(m+ 1) — (m + 1) lnp and 
U(a) = 2Z 2 aj{m + 1) — (m + l)ln<7. Ref. [12] shows that the superpartner 
equations and eigensolutions have the same form as Eqs. (6.1), (6.2) and their 
analogues for the a coordinate, but with the replacements n ] -4 n\ — 1, n 2 — + 


448 



n 2 — 1, m — *• m + 2. The supersymmetry leaves invariant the principal quantum 
number n, as needed for a mapping between degenerate states. 

The basic approach to a physical interpretation is analogous to the spherical- 
coordinate case but the existence of two supersymmetries permits more possibili- 
ties. Only two (natural) cases are mentioned here; see ref. [12] for more details. 

The valence-electron states for light alkali-metal atoms can be labeled using 
a spectroscopic notation based on the parabolic quantum numbers «i, n 2 , m. For 
convenience, the lowest- and highest-energy sublevels in a manifold of states with 
fixed m > 0 are called the red and blue levels, respectively. It turns out that 
a natural choice for the bosonic stack is the set of all blue states of hydrogen. 
Interpreting the absence of a fermionic-stack partner for the ground state as a 
consequence of the Pauli principle and treating electron-electron interactions as 
supersymmetry-breaking effects, the blue lines of hydrogen can be shown to be 
partnered by the blue lines of lithium. A similar construction can be made for 
the red lines of hydrogen and lithium. These supersymmetries extend to other 
alkali-metal atoms in a manner analogous to the spherical-coordinate case. 


7. Quantum-Defect Theory in Parabolic Coordinates 


The interactions breaking the supersymmetries in parabolic coordinates can 
be incorporated using notions from quantum-defect theory. For this case, the 
asymptotic quantum defects 8 entering the Rydberg expression for the eigenener- 
gies depend on the two parabolic quantum numbers n\ and m. It can be shown 
that an effective potential exists for the two-dimensional separated equation in p 
and a that has analytical eigenfunctions of the same form as in section 6 with the 
Rydberg series as eigenvalues, while incorporating the integer shifts characteristic 
of supersymmetry [12]. 

In separated form, the model is specified by 


r d 2 


1-mJ 

v 


p 2 


(7.1) 


and a similar expression for the equation in a with subscripts 1 replaced by sub- 
scripts 2. The functions x*(p) play a role analogous to Xi(p) in section 6, and the 
separation constants Z* and Z\ satisfy Z* + Z\ = 1. 

Solutions to these equations exist only if nj = Z* /y/—2E* — (m* + l)/2 and 
n 2 = /V—2E* ~ ( rn 2 + l)/ 2 are nonnegative integers. Note that n* = n — 8 = 

n i + n 2 + ( m i + m 2 )/2 + 1. The wavefunctions are 
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(7.2) 
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with an analogous result for the cr-dependent equation. It can be shown that these 
model wavefunctions are not fully orthogonal for physical values of the asymptotic 
quantum defect, although the off-diagonal matrix elements are numerically small. 

The four quantum numbers that appear depend on the parabolic quantum 
numbers, the quantum defect, and the supersymmetry integers. Define Uj = 
n i —lit = ^ 2 -/ 2 , where I\ and I 2 are integers, and raj = m+a 1 , = m + a' 2 , 

with the constraint c*i + 0-2 = 2 1\ + 21 2 — 28. If 8 is set to zero, the constraint 
reduces the four variables Jj, I 2 , aq, 02 to three and the exact supersymmetry 
limit of section 6 can be recovered. 


8. The Stark Effect 

In an external electric field of strength F, the energy levels of the valence 
electron in an alkali-metal atom are shifted by amounts depending on F. A plot 
of the eigenenergies as a function of F is called a Stark map. An exact theoretical 
treatment for Stark maps of alkali-metal atoms is impractical. However, numerical 
approximation methods yield Stark maps in good agreement with experiment [13]. 

The analytical eigenfunctions obtained from supersymmetry-based quantum- 
defect theories in parabolic and spherical coordinates provide another approach 
to obtaining Stark maps [12]. In parabolic coordinates, the nonorthogonality of 
the eigenfunctions is a barrier to exact analysis. One possible approximation 
is to neglect the small off-diagonal components of the zero-field energies. This 
calculation yields Stark maps with curves that are largely linear, in disagreement 
with experimental results. 

This difficulty is absent in spherical coordinates. The zero-field wavefunc- 
tions of the model described in section 5 can be used to calculate Stark matrix 
elements, and the ensuing energy matrix can be diagonalized for a subset of the 
perturbation basis. Ref. [12] contains detailed results for the n = 15 Stark states 
of lithium and sodium for m = 0 and m = 1 , evaluated over the ranges used in ref. 
[13] for the energies and electric field strengths. The non-hydrogenic structure is 
clearly reproduced, demonstrating the expected sizable anticrossings and small- 
field quadratic Stark effects for the s and p levels. All the Stark maps generated 
for the n = 15 lines of lithium and sodium by the supersymmetry-based theory 
are indistinguishable from the numerical and experimental results of ref. [13]. 
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